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, Abstract 

In this paper, we study the problem of multipath channel estimation for direct sequence spread spectrum signals. 
To resolve multipath components arriving within a short interval, we propose a new algorithm called the least- 
squares based iterative multipath super-resolution (LIMS). Compared to conventional super-resolution techniques, 
such as the multiple signal classification (MUSIC) and the estimation of signal parameters via rotation invariance 
techniques (ESPRIT), our algorithm has several appealing features. In particular, even in critical situations where 
the conventional super-resolution techniques are not very powerful due to limited data or the correlation between 
path coefficients, the LIMS algorithm can produce successful results. In addition, due to its iterative nature, the 

^-H \ LIMS algorithm is suitable for recursive multipath tracking, whereas the conventional super-resolution techniques 

may not be. Through numerical simulations, we show that the LIMS algorithm can resolve the first arrival path 
\ among closely arriving independently faded multipaths with a much lower mean square error than can conventional 

. early-late discriminator based techniques. 

Index Terms 

> ■ 

Multipath channel, super-resolution, pseudo-noise code, early-late discriminator, least- squares, maximum like- 
! lihood 

I. Introduction 

O ■ 

A pseudo-noise (PN) code sequence, which is a deterministic but noise-like sequence generated by a 
■ linear feedback shift register, is widely used in the areas of radar and sonar signal processing, ranging and 
a ^ ! positioning, and digital communication, since it is convenient to estimate the code phase of the received 
^ \ signal using its peaky auto-correlation function (ACF). When the PN code sequence is launched to a 
channel and received by a remote terminal, the cross-correlation of the received signal and a receiver 
replica of the PN code sequence yields a sharp peak at the lag of the propagation delay caused by the 
channel. This delay is related to the physical separation between the transmitter and the receiver, and 
thus gives important information for source localization or target detection. However, in many practical 
situations, the propagation medium suffers from multiple echoes with short delays due to scattering objects 
around the receiver. In such a case, we are frequently interested in estimating the delays of those echoes 
and, especially, the first arrival path, since it is the most likely to be the line-of-sight path, which contains 
the information of the true distance between the transmitter and the receiver. 

In the multipath channel, the conventional PN code phase estimation method using the ACF, which 
finds the peak of the correlator output, would probably give an incorrect estimate of the first arrival 
path delay unless the first arrival path is the most dominant path. This is because the ACF for the first 
arrival path is buried under the sum of other overlapping ACFs of echos, and thus the lag of the largest 
magnitude correlator output is often different from the first arrival path delay. We can work around this 
problem to some extent by finding multiple local extrema of the correlator output and selecting the one 
with the smallest lag as the estimate of the first arrival path [1J. However, the time resolution of this 
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method is limited to the chip duration of the PN code sequence, which is not sufficient in some narrow 
band applications. 

In order to improve resolution of the conventional ACF based code phase estimation, super-resolution 
techniques can be taken into account. The super-resolution techniques, such as the multiple signal classi- 
fication (MUSIC) J2]] and the estimation of signal parameters via rotation invariance techniques (ESPRIT) 
0, were originally proposed for direction finding of multiple targets with passive sensor arrays, and 
first applied to the channel multipath resolution in [4j. In H] and Q, the super-resolution techniques 
were applied directly to the received signal samples without cross-correlation. In particular, in [5], it was 
shown that the super-resolution techniques can be applied to received signals sampled at a considerably 
lower rate than the chip rate of the PN code sequence when the signals have a finite rate of innovation 
(FRI). On the other hand, in J6), Q, super-resolution techniques for the correlator output signal were 
introduced. These techniques can significantly improve the resolution of multipath delays, compared to 
conventional techniques, but have some major drawbacks that render them impractical. First, since the 
super-resolution techniques are based on the subspace decomposition of the signal correlation matrix, 
the signal components contributed by different channel paths should occupy distinct dimensions of the 
subspace in order to be distinguishable. For this, each distinct path should be far enough apart from 
the others, and the attenuation of each path should be varying randomly and not be perfectly correlated 
with the others during an observation interval. These conditions are not satisfied when the paths are quite 
narrowly spaced, the channel is slowly varying, or the observation interval is not long enough. In addition, 
it is not so reasonable to assume that only the channel attenuations are randomly varying while the path 
delays are fixed. As a result, additional techniques, such as frequency smoothing fl6), 0, are required at 
the transmitter to make the assumption valid. As a second drawback, most super-resolution techniques are 
batch algorithms that compute a set of estimates from a large number of stationary observations. Although 
the technique for FRI signals [5] can work with a very small number of observations, it has very low noise 
immunity and also needs a large number of observations to achieve robustness against noise. The radio 
ranging application, however, often requires tracking of the signal as well. That is, a new observation, 
which can be non- stationary, is periodically given, and a new set of estimates should be immediately 
computed from the new observation and the previous set of estimates in a recursive manner J5), B3- Thus, 
the super-resolution techniques are not adequate for applications where tracking is a concern. 

In this paper, we consider a least-squares (LS) approach to resolve the short-delay multipath components 
with a correlator output signal. When the noise is Gaussian distributed, our approach includes the maximum 
likelihood (ML) approach, which is optimal in performance. The LS and ML approaches for the channel 
multipath resolution were presented directly and indirectly in earlier studies [fTOl . [fTTTl . |[T2'1, but they have 
not gathered much interest since the highly nonlinear nature in the parameters to be estimated renders 
the use of these approaches computationally intensive. In ifTOll and [fTTTl . faster algorithms for solving a 
nonlinear LS problem were derived by exploiting the linear predictability of sinusoidal signals, but they 
are not applicable to our case. In our work, we propose an iterative algorithm to solve the LS multipath 
resolution problem, which we will call the least-squares based iterative multipath super-resolution (LIMS) 
algorithm. Thanks to its iterative and recursive mode of operation, it is appropriate for multipath tracking 
with sequential observations of the received signal. Moreover, we show that the tracking of a PN code 
phase using the early-late (EL) discriminator fl8), ||9), [fT3ll is a special case of our algorithm. Through 
numerical simulations, we show that our scheme asymptotically achieves the optimal performance given 
by the Cramer-Rao bound (CRB) at high carrier-to-noise density ratio. In addition, we show that our 
algorithm outperforms the EL discriminator based technique in the presence of severe multipaths and 
noise. 

This paper is organized as follows. In Section [III we introduce the signal model and the LS formulation 
for the multipath resolution. In Section [Till we present the LIMS algorithm as a technique to solve the LS 
multipath resolution problem, and discuss the application of the LIMS algorithm for multipath tracking. 
The complexity of the LIMS algorithm and its relationship with the conventional EL discriminator based 
technique are also discussed in Section [III Results of numerical simulation in various channel models 
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and a comparison with the conventional method are given in Section JV] and Section [V] concludes the 
paper. 

We will use the following notation throughout this paper. Vectors or matrices are denoted by boldface 
symbols. The p-th element of a vector and the (p, q)-th element of a matrix are denoted by [-] p and ['} p , q , 
respectively. The superscripts T and H denote the transpose and conjugate transpose, respectively. The 
Euclidean norm of a vector is denoted by || • ||, and the infinity norm of a vector, i.e., the largest absolute 
value of elements, is denoted by || - ||oo- The fields of real and complex numbers are indicated by R and 
C, respectively, and the TV x TV identity matrix is denoted by I N . The statistical expectation is denoted 
by E{-} and the real part of a complex value is denoted by ${{•}. Finally, O(-) is the big O notation. 



II. System Model 

Consider a multipath channel between a transmitter and a receiver with an impulse response 

K-1 

h{t) = ^27kS(t-r k ), 0<t<T„ (1) 

k=0 

where K is the number of paths, {70, . . . , 7k-i} are the unknown complex channel coefficients for each 
path, {t , . . . , tk-i} are the unknown arrival time delays for each path, S(t) is the Dirac delta function, 
and Ti (3> max{r , . . . , tk-i}) is a time interval of interest. We assume that the delays {r , . . . , tk-i} 
are distinct and r < T\ < • • • < tk-u without loss of generality. Thus, the delay of the first arrival path 
is t . The transmitter launches a continuous time PN code sequence s(t) with a chip duration T c to the 
channel (OQ), and the baseband received signal at the receiver is written as 

r)(t) = h(t) * s(t) + v(t) 

K-1 

= - r ^ + ^ 1 < r - (2) 

k=0 

where v(t) is a circularly symmetric, zero-mean, complex Gaussian noise process with auto-correlation 
N S((). For the first process at the receiver, the cross-correlation of rj(t) with s(t) is taken over the time 
interval and the correlator output is uniformly sampled at TV lags {0, T s , . . . , (N — 1)T S } to yield 

K-1 

Vn = yZ (nT s -r k )+w n) n = 0, . . . , N - 1, (3) 

k=0 

where R(() is the auto-correlation function (ACF) of the PN code sequence given by 

R(C) = ^ s(t)s*(t-()dt, (4) 

l i JO 

and w n is the filtered noise process w(() sampled at ( = nT s , where 



1 



Ti 



w(0 = 7F / v(t)s*(t-C)dt. (5) 



Ti Jo 

For a binary phase shift keying (BPSK) modulated PN code sequence with a sufficiently large integration 
time Ti > T C9 the ACF © is ideally given by 

jr + 1, ~T C < C < 0, 

R({) = {-4r + l, 0<C<T c , (6) 
otherwise. 
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For notational simplicity in the subsequent expressions, we denote y = [y , • • • , Hn-i\ t , cch = [70 > • • • 3 7k-i] t ? 
tcH = b~o, • • • j 7"x-i] T , and w = [w , . . . , ^at_i] t . Then, from (0), the covariance matrix of the noise 
vector, £" {ww^ }, is given by ^C, where C has elements 

[C\ M = ^-Eiw^w*^} = R((p - q)T s ), l<p,q<N. (7) 

From the above signal model, the LS estimation of the channel coefficients and the path delays is 
formulated as 

{c,t} = argmin||Gy- GA(t)c|| 2 (8a) 

c,t 

= argmin ||Gy - G (B(c, t)t + b(c, t))|| 2 , (8b) 

c,t 

where c = [c , . . . , c M -i] T E C Mxl and t = [t , ■ ■ ■ , t M -i] T E M Mxl are the vectors of channel 
coefficients and path delays to be estimated, respectively, and M is the assumed number of multipaths 
for the estimation. Since the true number of multipaths, K, is generally unknown, we assume that M can 
be different from K. Note that there are two equivalent LS formulations (l8ab and (l8blh and the related 
parameters are defined as 



[A(t)] re+1)m+1 = R (nT s - t m ) , (9a) 

Cm tm T c < ^ <^ tm 

% ¥ S S < n < (9b) 
0, otherwise, 

M-l 

b(c, t) = ^2 b ( c m, t m ), (9c) 
1 + , t - m ^ < n < ^, 

[b(c m , t m )] n+1 = { Cm .( 1 _i^y t t < n < ^ (9d ) 

otherwise, 

rc = 0, . . . , iV — 1, m = 0, ...,M-1, 

and G e C 7Vx7V is a proper nonsingular weighting matrix. For the least-squares estimation ([8]) to be 
feasible, we assume that TV > 2M, and A(t) and B(c, t) have full column rank. The first assumption 
can be reasonable by taking a sufficient number of samples, and so is the second assumption when the 
delays {£ , • • • , t M -i} are distinct and sufficiently spaced by a small fraction of T c lfT4l 2.4.7]. 

Remark 1: The signal model © is restricted to using the triangular shaped ACF © for simplicity, but 
it can be extended to general ACFs through a proper linearization around the sampling points. That is, 
for a general ACF R(C) 9 (l9bb and (l9cl) can be replaced by 

r-D/ +m dR(nT s -t m ) 
[b(c, t)] n+1 = £ c m ■ (i?(nT s - f m ) - dR ^s~t m ) \ (U) 

m=0 V ^ m / 

In addition, though we assume the uniform sampling with a sampling period T s in ©, all the arguments 
in this paper can be immediately extended to non-uniform sampling. 

Remark 2: If we let G = C~2, the weighting matrix G whitens the noise and thus the LS formulation 
® becomes the maximum likelihood (ML) formulation. In Section |IV] we investigate the performance 
with and without noise whitening by letting G = C"2 and G = I N , respectively. 
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III. Solving the Least-Squares Multipath Resolution Problem 

From the LS formulation d8j, we can observe that it is linear in c given t but nonlinear in t due to the 
dependency of B(c, t) and b(c, t) on t. Thus, to find the solution to the LS problem, a computationally 
expensive multidimensional search seems to be inevitable. However, there exists a useful workaround 
to avoid an exhaustive search: the gradient descent algorithm lfT5ll . The gradient descent algorithm can 
efficiently find the local minimum of the LS problem in an iterative manner, by stepping from the current 
guess of the LS solution to the direction that decreases the LS weight in ([8]) at every iteration. In the 
following subsection, we propose a gradient descent algorithm for the LS multipath resolution problem, 
which will be denoted by the least-squares based iterative multipath super-resolution (LIMS) algorithm. 

A. Least-squares based iterative multipath super-resolution algorithm 

For notational simplicity, let y = Gy, A(t) = GA(t), B(c,t) = GB(c,t), and b(c,t) = Gb(c,t) in 
®. The operation flow of the LIMS algorithm is described as follows. First, at the /-th iteration (/ is a 
positive integer), the intermediate guesses of the channel coefficient and the path delay vectors, cj_i and 
t/_i, are given from the previous (/ — l)-th iteration. Given t/_i, a new estimate of the channel coefficient 
vector, ci, is computed as 

d = (A H (t«_i) A(ti_!)) _1 A H (t,_!)y. (12) 

Note that the matrix inversion in (fT2l) exists with assumptions that G is nonsingular and A(t) has full 
column rank. The new estimate c/ (IT2l) minimizes the LS weight in (l8ab for a given tj_i, but it requires a 
matrix inversion that can sometimes be troublesome in terms of complexity or numerical stability issues. 
Thus, as an alternative, a refinement by the gradient descent algorithm can be used as 

d = cj_i + a • A^VO (y - A(ti_i)ci_i) , (13) 

where a is a properly chosen step size. The refinement c/ (IT3l) does not have complexity or stability 
problems like (fT2l) . but may converge more slowly than (IT2T) . 

For the second step, a refinement t/ is computed from t/_i and ci using the gradient descent algorithm, 
for which it is required to find the gradient of the LS weight in (l8bl) with respect to t. However, since 
the LS weight is nonlinear and non-analytic in t, it is not easy to compute the exact gradient. Thus, we 
compute an approximate gradient on the basis of the following assumptions: For a small perturbation 
vector e, 

B(c,t + e) ~ B(c,t), (14a) 
b(c,t + e) ~b(c,t). (14b) 

From (l9bl) , (l9ck and (l9dl) , the assumptions in (IT4l) are quite reasonable for a magnitude of e sufficiently 
small in the sense that HeH^ <C T s . Under these assumptions, it can be assumed that B(c/,t/_i) and 
b(c/,t/_i) are constant in a small region around t/_i. As a result, the gradient descent for t reduces to 

t, = t,_i + ft • 3? {B* (cj, t,_i) (y - (b(q, t|_i)t,_i + b(cj, t,_i)) ) } 

= t,_i + p • 3? [B H (a, t,_i) (y - A(Vx)q) } , (15) 

where ft is a proper step size. The step size ft should be small enough to satisfy the assumptions in (fT4l) . 
but a systematic criterion for ft that guarantees the convergence is not easy to find. 

Upon getting q using (IT2l) or ([T3T) . and ti using ([T5lh the iteration index / is increased by one and the 
same computations are repeated for c^ +i and t/ +i . In this way, the iteration continues until a specific stop 
criterion is met. There are a number of stop criteria for this kind of iterative algorithm, e.g., reaching 
a fixed number of iterations, thresholding the magnitudes of changes of variables, etc., and we omit a 
detailed discussion of these criteria. The LIMS algorithm is summarized in Table ID 
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B. Complexity issues 

The major burden of the LIMS algorithm is the matrix inversion in ([T2lh which is 0(M 3 ) of computation 
per iteration. The true number of multipaths, K, could be large in practice and, hence, the assumed number 
of multipaths, M, should be accordingly large. However, assuming that only a few of the multipaths are 
dominant, we can keep M small; say, M < 10. Therefore, when TV 3> M, the matrix inversion may 
not be too burdensome compared to the overall complexity. In addition, as noted in Section IIII-Al the 
matrix inversion can be avoided by using (fl~3T) instead, which is O(NM) of computation per iteration. The 
total number of iterations is another important factor for the complexity. However, as will be discussed in 
Section IIII-Cl the number of iterations can be considerably smaller when the LIMS algorithm is combined 
with a tracking function such as the delay locked loop. 

C. Multipath tracking 

In many practical realizations of radio ranging systems, a new observation of the received signal is 
obtained periodically. Successive observations may include different multipath channels when the channel 
is time varying. Under the assumption that the channels for the successive observations are not statistically 
independent, the time varying path delays can be efficiently tracked. One such tracking system is the delay 
locked loop (DLL) (S), (9). In the conventional DLL used for PN code phase tracking, the early-late (EL) 
discriminator is widely used to generate an error signal from the previous estimate of the path delay and 
the new observation of the received signal (correlator output). The error signal is then filtered by the loop 
filter in the tracking loop and a refined delay estimate is computed. Now, instead of the EL discriminator, 
we consider using the LIMS algorithm in the tracking loop. 

Suppose that the v-th observation of the received signal samples, is given. Using the LIMS 

algorithm, estimates of the channel coefficient vector and the path delay vector are computed. 
Then, when the (v + l)-th observation y^ +1 ) is given, the LIMS algorithm starts a new round of iterations 
with c^) and as initial values. As long as and are good estimates of the channel coefficients 
and the path delays, and the channels for the v-th and the (y + l)-th observations do not differ much, the 
new estimates c^ +1 ) and t^ +1 ) might be obtained in a small number of iterations. 

When the channel is not time-varying but corrupted by different realizations of noise at different 
observations, applying the LIMS algorithm with a small number of iterations for each observation can 
be seen as the application of the stochastic gradient method for the channel estimation lfT5l . Thus, as 
the number of observations increases, the accuracy of channel estimation can improve until it reaches the 
steady state. 



TABLE I 

Summary of the LIMS algorithm. 



Input 

y = [2/0, . . . iVn-i] t 



Known parameters 

Weighting matrix: G G C 
Step sizes: a and (3 



Initial values 
c and t 



Computation: I = 1, 2, . . . 




Definitions 

A(t), B(c,t): given in © 

y = Gy, A(t) GA(t), B(c,t) GB(c,t) 
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D. Relationship with the conventional early-late discriminator 

The EL discriminator computes the difference between the powers (magnitude squares) of two samples 
with different lags at the correlator output. The two samples are called the early and the late samples and 
are separated by a constant interval. When the peak of an ACF is located in between the two samples, 
the difference between the two powers serves as a feedback signal that shifts the sampling points so that 
the peak of the ACF is at the midpoint of the sampling points. Intuitively, this operation is established 
on the assumptions that the channel has a single path and that the correlator output maintains the exact 
triangular shape of the ACF R(Q (|6]). However, these assumptions are not valid when the channel has 
multiple paths and the received signal is corrupted by noise. 

Basically, the EL discriminator can be thought of as a special case of the LIMS algorithm with the 
single path assumption, i.e., M = 1, and two samples, i.e., TV = 2. To illustrate this relationship, let us 
first consider 2P sampling points at ( = , p = — P, . . . , P — 1, where P = 

sampling points, the sampled values of the correlator output are denoted by Ve(- p ) 
for p > 0. Likewise, the sampled values of the scaled ACF, c P(C — ^o), are denoted by r E (- p ) for p < 
and rc(p+i) for p > 0, as shown in Fig. [TJ Now, we temporarily let T s = T C9 and thus P = 1 and the 
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Note that, since the two samples are one chip apart, the noise levels for the two samples are independent 
according to ©. Hence, by Remark |2l (fT6l) is the ML weight as well. From ©, parameter vectors for 
the LS weight (fT6b are defined as 
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Fig. 2. Discriminator responses. 



Applying the LIMS algorithm to the LS weight (fT6l) with the initial value t ,o = 0, the first iteration gives 



c ,! = (A H (t 0) o)A(t ,o)) 1 A H (t 0fi ) 
by (O, and t ,i = t 0fi - (3D by (TBI) , where 



(18) 



D = -3?<|B if ( C o )1 ,to,o) 



Vli 



co,iA(t ,o) 



= -^{(Vei + Vli)* (vei ~ Vli)} 



Wei] 



\Vli\ 



(19) 



(20) 



As expected, D (fT9)) . the approximate gradient with respect to t , is the same as the EL discriminator. 
Similarly, for T s < T c , we can derive 

D = -SR {B H (co, u O)^ 1 (y - ^,^(0))} , 

where y = [y EP , . . . , y E1 , y LU . . . , y L p] T , 

c ,! = (A^OjC^AlO))- 1 A^(0)C-V, (21) 

and A(0) and B(c ,i,0) are defined similarly to (fT7i by letting M = 1 in ©. Note that, in (O and 
(I2TT) . noise whitening is applied. Then, (f20l) can be regarded as a generalized 2P-point discriminator 
as an extension of the conventional 2-point discriminator like the EL discriminator. The quantity D as 
a function of the code phase offset for different values of T s is evaluated and plotted in Fig. [2l The 
discriminator responses of the EL discriminators with sample spacing T c (wide discriminator) and 0.2T C 
(narrow discriminator) are also plotted in Fig. |2] for comparison. In the figure, it is observed that, as T s 
decreases, the discriminator curve gets steeper around the zero offset, which implies stronger resistance 
against perturbation after the DLL is locked at the peak of the ACF lfT3ll . 



IV. Numerical results 

This section presents simulation results to demonstrate the performance of the LIMS algorithm. For 
comparison, the conventional EL discriminators for PN code phase tracking with wide (T c ) and narrow 
(0.2T C ) spacings are also considered. In particular, performance is evaluated with the estimation accuracy 
of the delay of the first arrival path and, as the performance measure, the mean square error (MSE) of 
the delay estimate is used. For all simulations, the results are obtained by averaging over 10 4 trials with 



9 



independent realizations of noise and channel impulse response. As for the PN code sequence, we use 
the global positioning system (GPS) coarse-acquisition (C/A) code sequence lfT6j which has 1023 chips 
per period with a 1.023MHz chip rate; one chip occupies T c = 1/(1.023 x 10 6 ) seconds. The received 
signal is generated by superimposing the attenuated, phase rotated, and delayed replicas of the C/A code 
sequence according to the multipath channel realization, and adding white Gaussian noise with the power 
spectral density N . The received carrier power is defined as C = E {\h(t) * s(t)\ 2 } in (O, where the 
expectation is over the channel realizations. The carrier power C varies resulting in a carrier to noise 
density ratio (C/N ) between 20dBHz and 60dBHz. At the receiver, the integration time of the correlator, 
Ti, is set to be 10ms, which corresponds to 10 periods of the GPS C/A code sequence. 

For all simulations, it is assumed that a window of the correlator output signal of 3T C width is available 
for the EL discriminator and the LIMS algorithm. The center of the window is defined as the zero lag, i.e., 
( = 0. It is also assumed that the coarse acquisition of the PN code phase is perfect regardless of C/N 
so that the true first arrival path is placed at a random lag uniformly distributed over [— 0.5T C , 0.5T C ). 
Note that, since the size of the observation is too small to compute the sample correlation matrix of the 
received signal, conventional super-resolution techniques like MUSIC and ESPRIT cannot be applied to 
the considered simulation environment. 

A. Simulation system 

For the simulation of the LIMS algorithm, the sampling period T s = 0.1X C is used. As described in 
Section IIII-Al equation (fT2l) is used for the channel coefficient update and, for ease of implementation, 
the signed gradient descent algorithm is used in place of (031) for the delay update. That is, 



t, = t,_i + p • sgn [ft {b^(c«, t,_i) (y - A(t«_i)c«) } 



(22) 



where sgn[-] is the element- wise sign function. For the step size f3 in (l22l) . we use (5 = T s /64. To start 
the iteration, the initial value t = [to,o? • • • ?*m-i,o] t * s given such that i ,o == — 0.5T C , t M -i,o = 0.5T C , 
and £ m5 o, < m < M — 1 evenly divide the interval (— 0.5T C , 0.5T C ). With this initial value, the final 
estimates of the channel coefficient vector c = [c , . . . , c M -i] T and the delays vector t = [£ , • • • , i M -i] T 
are obtained by c 50 o and t 50 o, respectively. 

Upon completion of the estimation process, a simple path validation process follows to reduce the 
probability of false alarm; an estimated path is valid if the estimated power (|c m | 2 ) of the path is larger 
than 1% of the total estimated power (Ylm=l l^n| 2 )- The path with the smallest delay among the valid 
paths is chosen for the first arrival path. 

Finally, in the simulations of the LIMS algorithm, two different cases, without noise whitening, i.e., 
G = I M , and with noise whitening, i.e., G = C" 1 / 2 , are included for comparison. 



B. Two-path non-fading channel 

r i T 

The first test channel considered for simulation is a two-path non-fading channel with cch = 75? 75 
and tcH = [0, y] T - The performance results in terms of the MSE with respect to C/N are shown in Fig. 



3(a)[ As a benchmark, the CRB, which is derived in Appendix, is also plotted in the same figure. With 
this channel, the correlator output has a flat top when there is no noise. As a result, the conventional EL 
discriminator, which assumes a single path and the triangular ACF, cannot perform well with the channel. 
On the other hand, the LIMS algorithm can handle the two paths simultaneously and outperforms the 
EL discriminator. It is also observed that the LIMS algorithm with noise whitening consistently performs 
better than that without noise whitening. In particular, at high C/AT 's, the LIMS algorithm with noise 
whitening asymptotically achieves the CRB. 

In Fig. |3(b)[ learning curves of the LIMS algorithm with noise whitening for the first arrival path delay 
are depicted for 50 independent trials at C/N = 30dBHz. Observe that, even with a very large initial 
deviation, the algorithm successfully converges to the vicinity of the true first arrival path delay within 
300 iterations. 
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ON [dBHz] Number of iterations 

(a) Mean square errors (b) Learning curves at C/No = 30dBHz 

Fig. 3. Mean square errors and learning curves of the first arrival path delay estimation: cch = [l/y/2, 1/V2] and t C H = [0, T C /2] T . 




20 30 40 50 60 

C/N [dBHz] 



Fig. 4. Mean square errors of the first arrival path delay estimation: Single-path fading channel. 



C. Single-path fading channel 

For the second simulation, a channel with only one path is considered, where the channel coefficient 
70 is a zero-mean circular symmetric complex Gaussian random variable with i?{|7o| 2 } = 1- The MSE 
performances of the EL discriminators and the LIMS algorithms with different settings for this channel are 
plotted in Fig. HI As shown in the figure, the EL discriminators, especially the one with narrow spacing, 
perform very well for the channel. This is because the assumption that the correlator output has a triangular 
shape, on which the EL discriminator is based, is valid for this case. Among all performance curves plotted 
in Fig. HI the best is achieved by the LIMS algorithm with M = 1 and noise whitening. Compared to the 
performance of the EL discriminators, the best performance is an order of magnitude superior. In terms 
of the generalized discriminator presented in Section IIII-Dl this performance improvement arises from 
the fact that the LIMS algorithm can make use of a larger number of received signal samples and thus 
has a steeper discriminator response. Finally, in Fig. |4] note that the LIMS algorithm with M = 2 and 
noise whitening has a performance comparable to that with M = 1 and noise whitening even though the 
number of assumed paths is not the same as the true value. 
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Fig. 5. Mean square errors of the first arrival path delay estimation: Multipath fading channels without band limitation. 



D. Multipath fading channels 

The third test channel models considered for simulation are the multipath fading channels. Two multipath 
fading channel models, which are denoted by channels A and B, respectively, are considered. Channels A 
and B have three and four paths, respectively, and the power-delay profiles are given in Table HI Note that, 
for both channels, the average power of the first arrival path is 7dB lower than that of the strongest path. 
With the given power-delay profiles, the channel coefficients are generated independently as zero-mean 
circular symmetric complex Gaussian random variables with corresponding variances for each trial. The 
MSE performances of the EL discriminators and the LIMS algorithms with different settings are plotted 
in Figs. |5(a)| and |5(b)| for channels A and B, respectively. In Figs. |5(a)| and |5(b)[ it is observed that the 
EL discriminators yield very poor performances. This is because the EL discriminator tends to keep track 
of the largest magnitude point at the correlator output, which would be located mostly near the strongest 
path and hardly near the first arrival path. Due to the same reason, the performance of the LIMS algorithm 
with M = 1 is also bounded. However, the LIMS algorithms with M > 2 show much better performances 
since they can track multiple paths simultaneously, and one of the tracked paths is possibly close to the 
true first arrival path. 

The effect of limited bandwidth on the performance should also be investigated. In practice, the received 
signal is bandlimited to the pre-correlation bandwidth by a band pass filter before the cross-correlation 
with the PN code sequence replica, and thus the ACF may not have as sharp a triangular shape as @ 
has Q/71 4.2]. This causes a mismatch between the assumed signal model and the true signal, and can 
entail some performance loss. For the simulation, an ideal band pass filter with a rectangular frequency 
response centered at the carrier frequency of the signal is assumed. Figs. |6(a)| and |6(a)| show the results 
obtained for channels A and B, respectively, when an 8MHz pre-correlation bandwidth is applied. By 
comparing Figs. |5] and |6l it is found that there is no significant loss when using an 8MHz pre-correlation 



TABLE II 

Parameters of the three- and four-path fading channels. 



Path 


Channel A 


Channel B 


Relative power (dB) 


Delay (T c ) 


Relative power (dB) 


Delay (T c ) 


1 


-7.0 





-7.0 





2 





0.3 


-7.0 


0.2 


3 


-2.2 


0.5 





0.4 


4 






-2.2 


0.6 
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Wide EL (7) 
Narrow EL (0.27^ 

LIMS, w/o whitening, M=1 
LIMS, w/ whitening, /W=1 
LIMS, w/o whitening, M=2 
LIMS, w/ whitening, M=2 
LIMS, w/o whitening, M=3 
LIMS, w/ whitening, /W=3 
LIMS, w/o whitening, M=4 
LIMS, w/ whitening, M=4 



^_Wide EL (Tj 
9 - Narrow EL (0.2 TJ 
■H — LIMS, w/o whitening, M=1 
B - LIMS, w/ whitening, M=1 
■0— LIMS, w/o whitening, M=2 
$ - LIMS, w/ whitening, M=2 
^ — LIMS, w/o whitening, M=3 

LIMS, w/ whitening, M=3 
A — LIMS, w/o whitening, M=4 
A - LIMS, w/ whitening, M=4 




(a) Channel A 



(b) Channel B 



Fig. 7. Mean square errors of the first arrival path delay estimation: Multipath fading channels with 2MHz pre-correlation bandwidth. 



bandwidth. The results with a 2MHz pre-correlation bandwidth are also shown in Figs. |7(a)| and |7(b) 



There is apparent loss in this case, especially for the LIMS algorithms with M = 3 and M = 4 at high 
C/Nq's. However, the performance loss is still not significant, and, therefore, one can conclude that the 
LIMS algorithm is quite robust to the band limitation. 



V. Conclusion 

In this paper, the multipath channel estimation problem for direct sequence spread spectrum signals 
was investigated. A new technique based on the least-squares criterion, denoted as the least-squares 
based iterative multipath super-resolution (LIMS) algorithm, was proposed for the resolution of short 
delay multipaths. Unlike conventional super-resolution techniques such as MUSIC and ESPRIT, the LIMS 
algorithm does not rely on subspace decomposition, and, therefore, has practical advantages. In particular, 
the LIMS algorithm can work with a short observation of the received signal regardless of the correlation 
between path coefficients. In addition, due to its iterative operation, the LIMS algorithm is suitable for 
recursive multipath tracking. We also presented a view of the LIMS algorithm as a generalization of the 
conventional EL discriminator for PN code phase tracking. Through numerical simulations, it was shown 
that the performance of the LIMS algorithm is superior to that of the EL discriminator in estimating the 
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delay of the first arrival path. As a result, the LIMS algorithm can be a useful technique for applications 
such as radio ranging, global navigation satellite systems (GNSSs), radio propagation analysis, and time 
synchronization and multipath tracking in wireless communication systems. 

Appendix 

The Cramer-Rao bound for the first arrival path delay estimation 

From the signal model ©, the inverse Fisher information matrix [[T8] B.3.3] for the estimation of the 
channel coefficients and the path delays is computed as 

F _1 = (3? {D^C^D}) -1 , (23) 

where D e C Nx2K is a concatenated matrix defined as 

D = [A(tcH) B(c C H,tc H )], (24) 

and A(-) and B(-, •) are defined in ©. The mean square error (MSE) of each parameter estimation is 
lower bounded by the corresponding diagonal element of the inverse Fisher information matrix. Thus, the 
MSE of the first arrival path delay estimation is lower bounded by the (K + l)-th diagonal element of 
the inverse Fisher information matrix (1231) . 
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